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ABSTRACT 

Most tmciofluidic systems rely on one of two manners of 
fluid transport; pressure-drivBn or electrokinetically-drlveii 
flow. This investigation focuses on describing these flows in 
micro^ricated cbaimcls and small diameter capillary tubes. 
Plow characterization is accomplished by xnteirogation of 
tmcron-scale fluid regions through a powerful^ non-in<rusive 
flow imaging technique. Interesting phenomena have bee^ 
observed from these detailed examinations. Our results are 
presented in conjunction with an evaluation of mechanisms that 
potentially explain observed deviations from the Holmholtc- 
Smoluchowslci equation. In particular, wc show that observed 
perturbations of electrokinetio flow in open capilkrics might be 
caused by induced pressure gradients. We also show how 
ihese indxiced preesurc gradients may globally perturb the flow 
in an electroldiietically-driven microfluldic system. 



NOMENCIATURE 




Symbol 


Definition 


Units (SI) 


E 


electric field 


Vm^' 


B> 


axial electric field 


Vra' 


F 


Faraday constant (N^e) 


Cmol"^ 


Id 


Modified Bessel fiincdon 




L 


le&gdi of capillary 


m 


P 


Piceaauie 


Pa 


Q 


Volumetric flow race 


m s 


R 


oaiversal gas constant 


Jmor'K"* 


S 


source term 




T 


teoipcsrature 


K 


To 


wall temperature 


K 



U 


electroosmotic or 






electrophoretic vclodty 




Uo 


electrophoretiG velocity 
atTg 


m 8"' 


Co 


ion concentration 


mcl 


dV 


differential vohmic 




k 


thertnal conductivity 




r 


radial dimension 


m 


To 


capillary inner radius 


m 


u 


velocity distribution 




X 


axial dimension 


m 




axial position of change 
in wall zeta potcntiBl 


m 


y 


distance from sui&cc 


m 


z 


chargo number 






permittivity 






electric potential 


V 




Debye lengtli 


m 




viscosity 


Pas 


Po 


viBOOsi^ of water at To 


Fas 


Pb 


electric charge density 


Cm-' 


P 


resistivity 


obiam 




zcta potential 


V 




electroosmotic zeta potential 


V 


electrophoretic 2»ta potential 


V 



INTRODUCTION 

Originally emerging from the field of analyticat cbcmisfry^ 
mtcro total analysis systems OlTAS) have recently enjoyed 
wide interest. Micro total analysis systems and other 
miorofluidic systems lely on both pressure and electrokinetic 
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fxtechanisms for fluid transpcnt ConipliQatod relationships can 
arise between the physical characteristics of the mjcrochantieU 
and the behavior of the multi-component fluida flowing through 
the channelfi. These relaticnships have not been systcmalieally 
^died and are not yet completely Understood, fn practice, 
researcher? are forced to rely on the use of trial and error hi tbe 
design of miniaturized fLaidic systems. 

In this investigation, wc employed a caged fluorescence 
imaging techniciuo in an effort to dctcnnine the mechanisms 
that cODitribute to the fluid flow in microfluidic Systran^. Aided 
with characterizations of relatively simple flows, we endeavor 
to accurately predict the flow bdiavior in more conq}licatcd 
flows and geometries. 



BACKGROUND 

Electrokinetic transport refers to the combination of 
electroosmotic and electrophoretic transport For a historical 
review of electrokinetic Sieory, see Bursreeu and Nalcach 
(1964), Included below as a brief description of each type of 
transport. 



eiBctroosmosis 

Electroosmosis refers to ttiB bulk movement of an aqueous 
solution past a stationaiy solid sui&ce, due to an extemally 
l^tied electric field. Electroosmosis reqimcs fhc existence of 
6 charged double-layer at the soUd-liquid interface. This 
charged double layer results from an attraction between bound 
surface charges and ions in the passing fiuid. glass 
C^ltaries, surface silanoL groups become dcprotonatsd and^ 
therefbre, are negatively charged. This negatively charged 
surface OKracta positive ions present in flic flow. In Iha 
situations addressed in this paper, only a very thin layer near 
the wall has a net charge. 

Rice and Whitehead (1965} give a complete analysis of 
electroosmodo flow in roimd capillaries; here we quote only a 
few significant results to which we will reffcr later. At 
eqoilibrium, a Boltzmann distribution of the chaiges in the 
solution will CTcist (Probstein,- 1994) widi the cbaige density 
distribution represented by: 



In our experiments, ^i^ro is approximately 10*^; tirus tbe double 
layer is very OiiiL The equation of motion for steady, low 
Reynolds nombcr, incompressible flow is: 

;iV^u = -p£E+VP, (3) 

Id the limit of a small Debye length, X^, solving Eq. (3) yields 
the Helmholtt^oluchowki equation for the electroosmotic 
velocity. 

U^II£m1jl (4) 

The Helmholtz-Smoluchowski equation predicts a "plug- 
like" velocity profile when the Debye length is much less tiban 
the capillBiy diameter. The electroosmolic velocity^ U, is 
appioximat^y O.I nun/sec, when ^ = 0.1 V, 100 V/cm, 
uid fi is fhc viscosity of water (Probstein, 1 994), 



Electfophoresis 

Electrophoresis describes the motion of a charged sur&ce 
Submerged in a fluid under the action of an applied electric 
field. Consideting the case of a charged dyo molecule, it can be 
shown (Probstdn, 1994) that the electrophofetic velocity of the 
dye is again described by the Hehnholtz-Smolubhowski 
equation. 

U=^^^ (5) 

The electroosmotic zcta potential is a property of the capillary 
surface while the electroplioretic zeta potentisi is a property of 
the charged ctyo. In general, these two zeta potentials wiD not 
be equal. 

It is important to note that electroosmosis results in a net 
mass transfer of the aqueous solution; whereas, electruphoresis 
causes movement of charged particles or molecules trough a 
stationary soJution. For a netrueutral solution of charged 
molecules^ electrophoresis will not result m bulk motion of the 
solvent. Also note that both derivations shown above result m 
a velocity that is not a function of (he radial position; hence» it 
is termed a "plug-like" flow. 



The potential flcld is determined by solving the Poisson- 
Boltzmann equation, resulting in the following charge density. 



EXPERIMENTAL Dn^AlLS 

Two different types of microchanncls were used in this 
investigation. Initially, we used rectangular silicone 
mi^ochanncls. We then began investigating flows within 
round glass cq>illariefi with an internal diameter of lOQ |'Tn 
The capillaries were employed in an efibrt to isolate 
unexplained flow phenomena observed hi our initial silieone 
microcbaimel mvestigation. 

To create the silicone chanpels, silicon fmcroehannel 
molds were fabricated as ''negative" masters using standard 
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photoHthogcaphy. An STS deep reactive ion etch (DRIB) 
used to create raised rectangular stnictuies (invcitcd diannels) 
with dimeanozis of 10-200 in ividtb, 5 cm m length and 40 
in height An RTV silicone (Dow Coming Sylgard 1&4) 
Was used to form positive replicas of the silicon structures 
(Effenhauscr ct aL^ 1997). A photograph of a sealed replica is 
shown in Fig. ]. The elastomer channels were sealed wiih 
either glass microscope slides or elastomer coated ^lass 
microscope slides. 




Figure 1. Photograph of the elastomer device with three 
teet microchannols. Tho dreles are the Injection mile and 
the squares are used to study preeeure-drtven flow En a 
euddet) expansion. 

Our e^erimental approach utilized caged fluoresccnl dyes 
to image the flows. These dyes were activated (nncagod) by a 
UV laser beam and tracked by flooieecencc imagmg, see Fig. 2, 
(Paul et al^ 1998). In contrast to standard methods mjecting 
dye, this visualization technique allows deimitiaii of oarrow, 
fluorescent regions at virtually any location along the channel. 
Specifically, this allows for a precise definition of a small 
fluoresDont region and therefore an itnproved dctwminafion of 
the dye profile. Therefore, this caged fluorescence technique 
has advantages over injection methods that must transport ttie 
dye fiom the location of iDjection to the region of interest 

The caged dyes were purchased from Molecular Probes 
Inc. end used in millijnolar quantities. For the pressure-driven 
etndies, a caged fluorestscm was dissolved m distilled 




UV pulse initiates dye line 

I 

time evolved dye 

Figure 2. Top: Schematic of the experiment The Nd:YAG 
laser is used to uncage the dye and a Microblue 473 nm 
laser is used to excite the fluorescence. Bottom; 
Illustration of the uncaging process. The uncaging laser 
defines the starting fluorescent nuld volume. 

water. For the electroldnetic studies, a caged rhodamine was 
dissolved in Tris-EDTA buffer. The c^es are negatively 
cbazged when uncaged and may be cither chaiged or unchatged 
when caged. When caged, tho dye solution is noo-fluoresccnt 
Upon photo^tivation of a UV light illuminated volume, the 
protecting group is cleaved and the dye becomes fluorescent. 
Approximately 100 \iJ (per pulse) of 355 mn ligjit from a 
pulsed Nd:YAG laser nncaged the dye. The uncagod region 
extended roughly 20 along die length of the channel and 
spanned across its full widdi (Fig. 2). The laser beam was 
focused using a UV microscope objective that defined a shaip 
start zone. A continuous wave Microblue diode pumped laser at 
473 nm (Uniphase) was used for fluorescence excitation after 
activation (Fig, 2). Fluorescence images of the molecules 
dif&ising and moving in the local flow were collected using a 
microscope objective and a video-rate, interlaced camera 
(Texas Listiuments MC-780PH). 
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We have investigated transport in iziiciolluidi& systems and 
present expcrimeotal studies of pressuic-dnvcn and 
electroldnetiD flows. Wc also present numeric^ simulatloiu 
which were perfonned using the NctFloW module in Memcad 
(KIEMCAD, 1998). The NetFlow module employs a. diree- 
dimensional finite element based tool to solve iJic Navicr* 
Stokes equations. Wo note that tiie pressure-driven flow data 
show good agreement hctwcen (he preliminary experiments aud 
the simulations. Measurements of the elcctrokinetic flows, 
however, reveal deviations from simple tlieory. 

The evolution of a pressure-driven fluorescent profile is 
compBTod to a coircsponding BOqucnoe of simuladons In Fig. 3. 
The agreement is good, Imt note that there is dye retention near 
the walls that is not duplicated in the nmnerical simulations. 
This dye retention could be due to surface roughness. Atomic 
force microscopy indicates that the silicone channel has an 
RMS surfoce tougbness of 7 nm. 




Figure 3» Experimental (fefi) and simulated (right) pressuro- 
driven flow entering a sudden expanelon In an elastomer 
structure. The channel entering the expansion i« 
rectangular (100 \im wide and 40 deep) 

The evolution of an clectrokinedcally-driven fLuorescent 
profile in a round, fiised-silica capillary is shown in Fig. 4. 
These nnages reveal that the velocity profile appears to be 
constant across the channel, in agreement with £q. (4). Figure 
5 shows etectrokinetic flow in a DBwax-coatod (polyethylene 




Figure 4, Electrokinetic flow in a 75 (un diameter fiised- 
silica capillary exhibiting a plug-like velocity profile. 
Applied field strength was 350 volts/cm. The numbers with 
oach image corraspand to the time In mlUi$eoond9 after the 
uncaging pulse. 



Figure &. Electrokinetic flow In DBwax-coated (polyetfiyfene 
glycol), 100 pm diameter fkised-sllica capillary. The wax- 
coated capillaries were purchased from J&W Scientifio. 
Applied field strength was 245 volts/cm. The numbers with 
each Image correepond to the time In mllllseconda after the 
uncaging pulse. 

glycol) capillary. The wax coating is used to si^press 
electroosmodc flow, so that the movement of charged particles 
in the flow is due to electrophoresis only. However, these 
images show a pronounced parabolic velocity profile in the 
direction opposing the electric field. This profile does not 
agree with Eq. (5), and would be difficult to image without a 
high-resolution measuiemcnt technique such as caged 
florescence nnaging. Note that we observed no motion of the 
dye without application of the electric field, verij^ng that dieie 
was no externally applied pressure gradient Parabolic 
components to the electrokinetfc v^eeiiy profilo ha.vc 
previously been observed in capillaries (Paul et al., 1998) and 
inmicrochaonels (St. John et al., 1998). 
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DISCUSSION 

In Fig. 4, we see tliat the imaged velocity profiles for 
electfokinetic flow are in good agreement with Eq. (4). The 
flow is a combination of electro osmosis and electrophoresis, 
and the theory outlined above shows that both components 
yield a flat profile. Figure 5, however, shows significant 
deviationfl from a flat profile. This flow suggests a velocity 
pTofile similar to that observed in a ptessure-driven flow. 

In the following sections, we explore two mechanisms that 
potentially cxplabi observed departme from phig-like flow: 
vaiiatioiis in viscosity due (o temperature gradients, and 
Induced prcflsure gradients caused by a non-unifbrni wall 
potential (jLe. non-imifbim wax coating in the ease of Pig, 5)- 



The predicted velocity profiles for both an open capillaiy 
and a DBwax-costcd capillary each show an expected radial 
velocity variation of less than 0,1% (Pig. 6). The radial 
vcloci^ distdbutioa is parabolic in nature, but its magnitude is 
too small to explain the dye proflles seen in Fig, 5, To obtain a 
signxflcant radial variation in the electrophoretic velocity, one 
would need to increase the field strength or the buflbr 
concentration considerably. Our temperature variation model 
and, hence these conchisions, are in closd agreement imtfa those 
presented by OiwhJka et al, (1989), 



VarlatfonofvlficoeHv 

In this section, we examine the possibility that heat 
generation wiliiin the capillaiy could cause the electrophoretic 
veloGity profile to deviate from plug-like. To this end» we 
investigate the temperature^ viscosity and corresponding 
velocity distribudons across die capillary. 

In cylindrical coordinatesj we write the tenqieratUTe 
governing equation for our system as: 

rdr or p 

The source term, S, has units of power per unit volume. In this 
analysiSp we assume p Tcmains a constarit (or at least is not a 
strong function of temperature) based on the cjqierimentally 
determined fluid-filled capillary resistance and oapfllaiy 
dimensions. This empiiically determined resistivity agrees 
closely with those tabulated in the literature (Masiiyah 1994). 
Integrating twice and applying the boundary conditions of finite 
ten^eratute at the ccntcrline and constant wall temperature^ 
we obtain the tempezature distrTnitjon 

As an bitial approximation, we assume that the fluid within our 
capillary has the same thennal propertied and viscosity as 
water, As such, we have chosen to model the viscosity as a 
function of temperature, by: 

^=«q>[a+i(^)+c(^)'] (8) 

with T« " 293.15K and the constants given by a = -1.94, 
b - *4.80j and c « 6.74, as suggested by Reid (1 987), 

In order to relate the ten^eratureHdependent viscosity Co 
the electrophoreitio velocity profile withhi (he cftpillazy, we 
tetum to orUT analysis of electrophoresis. Utilizing Eq, (5) and 
negjecting the tenqyerature dependence of e, ^, or E,, we airive 
at 




Figure 9- Efoctrophoretlc velocity distribution for an open 
capillary and a DBwax-coated capillary 

From this analysis, wc conclude that the viscosity and, 
therefore, vdocity do not vary significantly across the capillary 
for the prescribed small-bore capillary parameters, buffer 
concentrations, and field strengths. Theicfrirei Joule heating 
induced veloeily variation does not appear to be a feasible 
mechanism frir fiie observed deviations fix>m plug-like flow. 



Induced Pressure Gradients 

In deriving Eq-(4), the pressure gradient term was 
neglected because there were no extartalfy applied pressure 
gradients. However, this expression will not bo correct if the 
electrokinetic transport induces a pressure gradicnL We now 
explore how the electroktoetic flow could potentially induce 
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pressure ^dicnts. If wc take the divGrgcnce of Eq. (3)j and 
apply contimiity fbr the fluid, (V tt)=o, we dBrive the 
foliowing expression. 

V^i> = E-Vp^+(V-EV^ (10) 

Equation (10) demonstrates how the pressure field evolves 
from gradientB in the charge density and applied field. 
Referring to Bq. (2)^ w& see that the fiiBt tenn on the ri^t hand 
aide of £q. (10) will be non-zero if <hc wall potmtial changes in 
Che axial direction. Nonuni&raiitics in the wall potential could 
exist due to eontatnination in the capiUary^ variations in wall 
eoatingSj or gradients in the buffer pH. Using classical 
techniqaBs for solving Poisson's equation (e.g. Green's 
fimcdons), it is possible to develop a general solution to Eq. 
(10) for an arbitiaiy zcta potential distribution along the wall of 
a cylindticaL capillary. Tn this paper, a direct method will be 
tised Co solve &r the pressure field in a capiUaty with a wall 
potondal as shown in Fig. 7. 



Region. 1 Regigu 2 




Figure 7. 3chennatlc of capillary with stop change In waH 
potenUal. The atiKcipatod velocity profiles are abo shown. 



capillary is a known reference value (i.e. atmospheric pressure). 
For an incompressible liquid, the volumetric flow rates must be 
the same in regions I and 2. Furthermore, we know that the 
pressure is continuous and piecewisc linear since wc assumed 
that the velocity prolile doc£i not vary axially in cither region 
(i.e. the pressure gradient must be constant in both regions). 
Therefore, we arrive at the following expression fbr the 
pressure in the capillaiy, where x is the axial position and Xq is 
the axial location of the discontinuity in wall potential. 




X 



Figure 8. Non-dimnnefonal pressure in a capillary with a 
step change In wall potential 



Let US assume that the fhiid velocity is only a function of 
radial position in both Region 1 and Region 2 shown in Fig. 7. 
Nate that Region 1 describes a section of the capillary with zero 
wall potential; whereas^ Region 2 has a finite wall potentiaL In 
other words, we will assume that axial velocity gradients are 
confined to a vanishingly fimall region near the discontinuity in 
wall potentiaL \^iiih this assumption, the volumetric flow rates 
in regions one and two are shown below. 



(11) 



We eenalder the ease whete the pctssuiB at eHbsc end of the 



(12) 



The presauro field corresponding to Fig. 7 is plotted in Fig. 
8. We see that the pressure gradient is constant in each section 
and acts to ensure that the mass flow rates are equal in each 
rcgioiL Figure 9 shows the predicted velocity profiles in 
regions 1 and 2. The velocity in the first section is due to 
pressure alone, since the wall potential there is zero. The 
velocity profde in the second section is simply a superposition 
of the flat electroosmotic profile and (he Poiseuille profile due 
to the iiiduccd pressure gradient The prespure gradient in 
each region is readily calculated from Eq. (12). The velocity 
profiles in Fig. 9 were also shown schematically in Fig. 7. 
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Figure 9. Radial velocity profiles in region 1 and region 2. 

To understand the sohition qijalitatively, one can imagine 
that the electroosmotic flow in Region 2 "pulls" the flow in 
Ff^on 1. In other wordSp the elecn-oosmotic nxovement of 
liquid in Region 2 tendfi to create svetion « the interface 
hetween the two regions. Hence, the presswe drops at tlie 
inter&co until the mass flow rates in each region ate equal. 
Note that the pressure-driven component submicis fixsm die 
electroosmotic flow in Region 2 whUe only pressure-driven 
flow is present in Region 1. The analytically calculated 
velocity profiles also agree qualitatively with the numerical 
simulation shown in Fig. 10 (MBMCAD, 1998). The 
simulation confiims thflt axial gtudicnts m the velocity field arc 
confined to a region less than a single diameter from the wall 
potential di$C0Dtiaiiity, 

Supeipositjons of Eq. (12) can be vised to find the pressure 
field caused by one or many small regions of 2ero wall 
potential. This is then a model of a "dirty" silica ceqiillary. 
ngnre 11 demonstrates the nondimcnslonal pressure field in a 
C^Uaiy with three ''dii^' regions, each covering 1% of the 
oapiUaiy's length. 

Note that flie pressure gradients in the clean regions are 
identical and that fhe locaKzcd di^ontimiitics in wall potential 
perturb the flaw everywhere. Equation (12) can also be 
snpczposcd to find die pressure fidd in a coated capillary (zero 
wi^ potential) with sm^ ^lean" regions where the wait 



Region 1 R^on 2 
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Figure 10. Simulation of Row in a capillary with a elep 
change In wall potential. 
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Figure 11. NondlmenGional pressure field in a capillary 
with three ■'dirty" spots 

potential is pon-zeco. Figure 12 shows the resulting velocity 
profile in the coated regionSj for various amounts of exposed 
(vncoat^) capillaty. The induced pressure gradient model 
presents a rational explanation of the fiow imaged in Fig. 5, as 
it admits a parabolic profile for reasonable variations in wall 
potential. 
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1 exposed re^a(19i rooJ) 
10 caiposed rcglolu (10S4 loui) 



Figure 12. Velocity profile in a coated rogion of a wax 
coated capillary vnth small regions uncoatad. The wax 
Coating acts to prevont elQctroo&niotJc flow. 

CONCLUSIONS 

We have imaged pressure-driven and electrokinetically- 
drivcn flows in silicone microciiannels and sitica capillaries 
using a powerful caged fluorescence technique. The pressuic- 
diivcn flows agree witli theory and numerical simulations. 
Single theory predicts that the eleclrokinetic flows should 
display unifonn (plug-like) velocity profiles. However, in 
some caseSk we ohsenre latge parabolic-likc components to the 
electrokinetic velocity profiles. 

We explored two meehanisnis m an attempt to explain the 
observed velocity profiles: viscosity variations and induced 
pressure gradients. We have demon^ated that Joule heating 
find the subsequent temperature and viscosity gradients do not 
adequately explain the observed parabolic veloci«y profile, but 
that induced pressure gradients are a possible explanation. An 
important consequence of the induced pressure gradient 
analysis is th&t small regions of non-uniform y^^all potential 
create pressure gradients everywhere In (he Jlow. This 
conclusion has important implications for those wishing to 
build lab-on-a-chip devices, since pressure gradients tend to 
mcrease the dispersion of sample packets in the flow. We are 
cuzreatly per&miing experiments that investigate the proposed 
induced pressure gradient mechatusm. 
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